November 28, 2012 3:16 WSPC/INSTRUCTION FILE SUSY 



International Journal of Modern Physics A 
© World Scientific Publishing Company 



Theoretical Overview on the Flavor Issues of Massive Neutrinos 



SHU LUO 

Department of Physics and Institute of Theoretical Physics and Astrophysics, 
Xiamen University, Xiamen, Fujian, 361005 China 
luoshu@xmu. edu. cn 



ZHI-ZHONG XING 

Institute of High Energy Physics and Theoretical Physics Center for Science Facilities, 
Chinese Academy of Sciences, Beijing 100049, China; 
Center for High Energy Physics, Peking University, Beijing 100080, China 
xingzz@ihep.ac. cn 



We present an overview on some basic properties of massive neutrinos and focus on 
their flavor issues, including the mass spectrum, flavor mixing pattern and CP violation. 
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a TeV seesaw scenario. 
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1. INTRODUCTION 

It is well known that two important experimental results were in the news in 2012: 

• On 8 March 2012, the Daya Bay Collaboration announced a 5.2a discovery 
of 9 13 =/= for this smallest neutrino mixing angle 

sin 2 20 13 = 0.092 ± 0.016(stat) ± 0.005(syst) (±lcr range) , (1) 

which is equivalent to 13 ~ 8.8° ± 0.8°. The convincing Daya Bay result 
puts the preliminary results of T2K ^ MINOS ^ and Double Chooz S 
which all hinted at 9 13 ^ in 2011, on solid ground. In particular, the 
fact that 8 13 is not strongly suppressed is a good news to the experimental 
attempts towards a measurement of CP violation in the lepton sector. 

• On 4 July 2012, the ATLAS 5 and CMS © Collaborations at the Large 
Hadron Collider (LHC) independently announced the discovery of a Higgs- 
like boson at the mass scale of 125 GeV to 127 GeV. If this result turns out 
to be true, it will have an important impact on the development of neutrino 
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physics because most of the neutrino mass models depend on the existence 
of the Higgs particle(s) and Yukawa interactions. 

Therefore, a brief overview of where we are standing and where we are expecting 
to go makes sense. 

The remaining parts of this review paper are organized as follows. In section 2 
we give a fast overview of some fundamental neutrino properties, such as the speed 
of neutrinos, the nature of massive neutrinos and the number of neutrino species. 
Section 3 is devoted to a brief description of the flavor issues of charged leptons and 
neutrinos, including the mass spectrum, flavor mixing pattern and CP violation. 
We compare the observed pattern of quark flavor mixing with that of lepton flavor 
mixing. In section 4 we go into details of possible lepton flavor structures by outlining 
two phenomenological strategies and taking a number of typical examples. The 
impact of large 9 13 on the running behaviors of other flavor mixing parameters 
is discussed in section 5 by using the one-loop renormalization-group equations 
(RGEs) in the framework of the minimal supersymmetric standard model (MSSM) . 
Section 6 is devoted to the seesaw-enhanced electromagnetic dipole moments of 
three Majorana neutrinos based on a TeV seesaw scenario. A summary and some 
concluding remarks are given in section 7. 



2. IMMEDIATE QUESTIONS ON NEUTRINOS 
2.1. Really Superluminal? 

The constancy of the speed of light c in vacuum and the independence of physical 
laws from the choice of inertial systems are two fundamental propositions of the 
special relativity (SR) ^. If our world is Lorentz invariant, a free particle's energy 
E, momentum p and rest mass m satisfy the relationship y^E 2 — \p\ 2 c 2 — mc 2 . 
The velocity of this particle turns out to be v = c^/l — m 2 c^/E 2 , implying that it 
cannot travel faster than light in vacuum. Could a particle be superluminal? The 
answer would be yes if the particle had an imaginary mass (called a "tachyon" ®) 
or if the Lorentz invariance were broken. 

The OPERA Collaboration claimed a "convincing" measurement of the super- 
luminal neutrinos in September 2011^. But five months later this story ended up 
with a mistake of the bad connection of the optical fiber. The OPERA paper was 
updated in July 2012 by including the new sources of errors, and the new result 
was in agreement with the SR. Here let us quote Steven Weinberg's comments on 
the original result of the OPERA experiment: "The report of this experiment is 
pretty impressive, but it bothers me that there is plenty of evidence that all sorts 
of other particles never travel faster than light, while observations of neutrinos are 
exceptionally difficult. It is as if someone said that there are fairies in the bottom 
of their garden, but they can only be seen on dark, foggy nights." 

An early measurement of the neutrino speed was done by using the pulsed pion 
beams (produced by the pulsed proton beams hitting a target) at the Fermilab in 
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Fig. 1. A schematic illustration of the "flavor hierarchy" and "flavor desert" problems in the 
fermion mass spectrum at the electroweak scale. Here the masses of three neutrinos are assumed 
to have a normal hierarchy. 



the 1970s UnUll l n this experiment the speed of muons was compared with that of 
neutrinos and antineutrinos. The same measurement was repeated in 2007 by using 
the MINOS detector In 2011 the speed of neutrinos was also measured in a 
few other lon g-ba seline neutrino experiments, such as the ICARUS tt3HH ; Borexino 
0^ and LVD experiments. But the most stringent constraint o n the spe ed of 
neutrinos was from the observational data of the Supernova 1987A H 7 | 18 | 19 [ | v _ 
c\/c < 10~ 9 obtained by comparing the arrival time of light with that of neutrinos. 



2.2. Definitely Massive? 

The neutrinos are massless in the standard model (SM) as a result of its simple 
structure and renormalizability. On the one hand, the SM does not contain any right- 
handed neutrinos, and thus there is no way to write out the Dirac neutrino mass 
term. On the other hand, the SM conserves the SU(2) L gauge symmetry and only 
contains the Higgs doublet, and thus the Majorana mass term is forbidden. Although 
the SM accidently possesses the (B — L) symmetry and naturally allows neutrinos 
to be massless, the vanishing of neutrino masses in the SM is not guaranteed by any 
fundamental symmetry or conservation law. Today we have achieved a lot of robust 
evidence for neutrino oscillations from solar, atmospheric, reactor and accelerator 
neutrino experiments. The phenomenon of neutrino oscillations implies that at least 
two of the three neutrinos must be massive and the lepton flavors must be mixed. 
This is the first convincing evidence for new physics beyond the SM. 

Fig. 1 is a schematic plot of the mass spectrum of the SM leptons and quarks 
at the electroweak scale. One can see that the span between m 1 and m t is at least 
twelve orders of magnitude. Furthermore, there exists an obvious "desert" spanning 
six orders of magnitude between the neutrino masses and the masses of the charged 
fermions. Why do the SM fermions have such hierarchy and desert puzzles? The 
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answer to this important question remains open. In particular, t he tiny neutrin o 
masses must have a peculiar origin (e.g., via the seesaw mechanisms ^0 | 21 | 22 | 23|24 fy 
Moreover, there might exist one or mor e keV sterile neutr inos in the desert as a 
natural candidate for warm dark matter 125126 27 28 29 30 SI] 



2.3. Dirac or Majorana? 

A pure Dirac mass term added into the SM is in general disfavored, unless the theory 
is built by introducing extra dimensions. Such a mass term in a renormalizable model 
of electroweak interactions would worsen the fermion mass hierarchy problem. An 
effective Majorana mass term given by the right-handed neutrinos and their charge- 
conjugated counterparts is not forbidden by the SM gauge symmetry, unless the 
contrived assumption of lepton number conservation is imposed on the theory. Hence 
most theorists believe that massive neutrinos are more likely to be the Majorana 
particles and their salient feature is lepton number violation. 

The unique window to verify the Majorana nature of massive neutrinos is to 
observe the neutrinoless doube-beta (0ufij3) decay. So far we have not obtained 
very convincing evidence for this lepton-number- violating process. Even if the 0^/3/3 
decay were never observed, one woul d still be unable to conclude that massive 
neutrinos are the Dirac particles I 32 | 33 [ -j^g effective mass of the Ov/3/3 decay could 
vanish if the Majorana CP-violating phases lie in some specific regions. On the other 
hand, there are some other mechanisms which can lead to the 0v(3f3 decay. Such 
new physics effects could be o f the s ame order as or even larger than the standard 
light-neutrino-exchange effect I 34 | 35 [ 

Given the SM interactions, a massive Dirac neutrino can have a tiny (one-loop) 
magnetic dipole moment /i„ ~ 3 x 10~ 20 /i B (m„/0.1 eV), where /j, b is the Bohr 
magneton I 36 | 37[ j n con trast, a massive Majorana neutrino cannot have magnetic 
and electric dipole moments, because its antiparticle is just itself. Both Dirac and 
Maj ora na neutrinos can have transition dipole moments (of a size comparable with 
P"v) which may give rise to neutrino decays, scattering effects with electrons, 
interactions with external magnetic fields (red-giant stars, the sun, supernovae, and 
so on), and contributions to neutrino masses. Current experimental bounds on the 
neutrino dipole moments are at the level of \i v < a few x 10 _11 /i B . 



2.4. More than Three Species? 

It is well known that "three" is a mystically popular number in particle physics, 
such as three Q = +2/3 quarks, three Q = —1/3 quarks, three Q = — 1 leptons, 
three Q = neutrinos, three colors and three forces in the SM. In this case, why do 
people want to go beyond N v = 3? 

The stu dy of light sterile neutrinos has become a popular direction in neutrino 
physics ^21. One is motivated to consider such "exotic" p articles for sev eral rea- 
sons. On the theoretical side, the type-I seesaw mechanismJ ^ 1 2 1 1 2 2 1 2 3 1 24 ] p rov ij es 
a very elegant interpretation of the small masses of v i (for i = 1,2,3) with the 
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help of two or three heavy sterile neutrinos, and the latter can even help account 
for the obs erve d matter-antimatter asymmetry of the Un iverse via the lep togenesis 
mechanism SHI. On the experimental side, the LSNdSH, MiniBooNE S21 and reac- 
tor E3I antineutrino anomalies can all be explained as the active-sterile antineutrino 
oscillations in the ass umptio n of one or two species of sterile antineutrinos whose 
masses are below 1 eV^ ' '. Furthermore, a careful analysis of the existing data on 
the Big Bang nucleosynthesis HH1 or the cosm ic microwave background anisotropy, 
galaxy clustering and supernovae la I 47 | 48 | 49 l seems to favor at least one species of 
sterile neutrinos at the sub-eV mass scale. On the other hand, sufficiently long-lived 
sterile neutrinos in the keV mass range might serve fo r a good candidate for warm 
dark matter if they were present in the early Universe ^1 

If the three known neutrinos have mixing with a few new degrees of freedom 
above or far above the Fermi scale, an exciting window will be open to new physics 
at high energy scales. In this case, however, the mixing between light and heavy 
neutrinos violates the unitarity of the 3x3 light neutrino mixing matrix and might 
result in some observable effects in the future precision neutrino experiments ^4 

3. NEUTRINO MASSES AND FLAVOR MIXING 

There are three central concepts in flavor physics: mass, flavor mixing and CP 
violation 1^21. The phenomenon of lepton flavor mixing at low energies is effectively 
described by a 3 x 3 matrix U, the so-called Maki-Nakagawa-Sakata-Pontecorvo 
(MNSP) matrix E2 | 53 | Qi ven the unitarity of U, it can be parametrized in terms of 
three angles and three phases 1511 



^ — S 12 C 23 c 12 S 13 S 23 el15 C 12 C 23 s 12 S 13 S 23 eZl5 C 13 S 23 1 (2) 




where = cosO^, = (for ij — 12, 13,23), and P u = Diag{e lp , e lCT , 1} is 

physically rel evan t if massive neutrinos are the Majorana particles. 

Fogli et aZ^Shave recently done a global analysis of current neutrino oscillation 
data and obtained the ranges of two neutrino mass-squared differences (dm 2 = 
m\ — m\ and Am 2 = |m§ — (m 2 + m 2 )/2|) and three neutrino mixing angles, as 
listed in Table 1, where NH and IH stand for the normal hierarchy (m 1 < m 2 < m 3 ) 
and the inverted hierarchy (rn 3 < m 1 < m 2 ), respectively. 

3.1. Neutrino Mass Spectrum 

The two mass-squared differences of three neutrinos have been determined, to a 
very good degree of accuracy, from current experimental data: Am|x = Sm 2 m 
7.5 x 10" 5 eV 2 and Am§ 2 w Am 2 w ±2.4 x 1CT 3 eV 2 . The absolute neutrino 
mass scale remains unknown and may hopefully be determined in the following 
experimental or observational ways: the single /3 decay, the Oz//3/3 decay and the 
cosmological constraints. Fig. 2 shows the parameter space of S = m 1 + m 2 + m 3 , 
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Tabic 1. Results of the global 3v oscillation analysis by Fogli et al in 2012, including the best-fit 
values and allowed la, 2a and 3a ranges for the neutrino oscillation parameters. 



Parameter 


Best fit 


la 


range 


2a range 


3a 


range 


<5m 2 /lCT b eV 2 (NH or IH) 


7.54 


7.32 


- 7.80 


7.15 - 8.00 


6.99 


- 8.18 


sin? 0i 2 /lO _1 (NH or IH) 


3.07 


2.91 


- 3.25 


2.75 - 3.42 


2.59 


- 3.59 


Am^/lO" 11 eV 2 (NH) 
Am 2 /10" 3 eV 2 (IH) 


2.43 
2.42 


2.33 
2.31 


- 2.49 

- 2.49 


2.27 - 2.55 
2.26 - 2.53 


2.19 
2.17 


- 2.62 

- 2.61 


sin 2 0i 3 /l(r 2 (NH) 
sin 2 6»i 3 /10- 2 (IH) 


2.41 
2.44 


2.16 
2.19 


- 2.66 

- 2.67 


1.93 - 2.90 

1.94 - 2.91 


1.69 
1.71 


- 3.13 

- 3.15 


sin* fes/lO- 1 (NH) 
si^Sas/lO" 1 (IH) 


3.86 
3.92 


3.65 
3.70 


- 4.10 

- 4.31 


3.48 - 4.48 
3.53 - 4.84 e 5.43 - 


3.31 
6.41 3.35 


- 6.37 

- 6.63 


8/tt (NH) 
S/n (IH) 


1.08 
1.09 


0.77 
0.83 


- 1.36 

- 1.47 









nip — \/m 2 \U el \ 2 + mWU e2 \ 2 + m 2 \U e3 \ 2 (the effective electron neutrino mass in 
the j3 decay) and = \m{U^y + m 2 U 2 2 + m 3U^ 3 \ (the effective mass of the Oi//3/3 
decay). A precision measurement of and £ in the sub-eV range could determine 
the neutrino mass hierarchy. In the two lower panels of Fig. 2 there remains a large 
vertical spread in the allowed slanted bands, as a result of the unknown Majorana 
CP-violating phases in the components. This observation indicates that more 
precise data in either the (jripp, Trip) plane or the (m^p, £) plane might provide 
some useful constraints on the Majorana phases. 

Before the absolute mass scale is determined, there remain two open questions: 
(1) is m 3 bigger or smaller than m-^l (2) can one neutrino mass (m 1 or m 3 ) be 
vanishing or vanishingly small? The first question awaits an experimental answer 
in the foreseeable future, such as a long-baseline neutrino oscillation experiment 
with appreciable terrestrial matter effects^] or a long-baseline reactor antineutrino 
oscillation experiment with accurate information on the energy spectrum | 57 | 58 | 
theoretical answer to the second question is strongly model-dependent. Examples of 
this type includ e the minimal type-I seesaw mechanism with two heavy Majorana 
neutrinos I 59 | 60 | or the Friedber g-Lee ansat z with an effective Dirac or Majorana 
neutrino mass operator IM162| 63 | 64 | 65 | 66 | 67 | 

3.2. Flavor Mixing Pattern 

The fact that the smallest neutrino mixing angle 13 is not strongly suppressed 
leads us to some new questions about the feature of lepton flavor mixing: (1) Can 
the relatively large 9 13 be understood by an underlying flavor symmetry or is it 
generated by a symmetry breaking mechanism or quantum corrections? (2) Does 
6*23 ^ 45° still hold? (3) What is the strength of leptonic CP violation? 

The structure of the MNSP lepton flavor mixing matrix U is significantly dif- 
ferent from that of the Cabibbo-Kobayashi-Maskawa (CKM) quark flavor mixing 
matrix V. The CKM matrix is nearly the unit matrix up to some small corrections, 
while the MNSP matrix has an approximate fi-T symmetry. The full fi-T symmetry 
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£ (eV) m p (eV) 

Fig. 2. Constraints obtained by Fogli et al in 2012 (at the 2o level) in the planes charted by 
any two among the absolute mass observables (the effective mass of the decay), rriap (the 
effective mass of the Ov/3/3 decay) and S (the sum of three neutrino masses). The blue (red) bands 
refer to the normal (inverted) neutrino mass hierarchy. 

of U in modulus is described by the equalities 

1^1 = 1^1, K 2 \ = \U T2 \, \U,a\ = \U T3 \, (3) 

equivalent to two independent sets of conditions in the standard parametrization 
given in Eq. (2) 

^23 = 45°, 13 =O°, (4) 

or 

6> 23 = 45° , S = ±90° . (5) 

If #23 is exactly equal to 45°, then one may arrive at a partial \x-t permutation 
symmetry in the MNSP matrix U (i.e., the equality \U„ 3 \ = \U t3 \). 

Now that 8 13 7^ 0° has firmly been established by the Daya Bay experiment ^ it 
becomes crucial to check the deviation of 6 23 from 45° and (or) a possible departure 
of (5 from ±90°. We speculate that U might have an approximate ji-r symmetry with 
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\U j\ ~ \U Ti \, in contrast with the approximate off-diagonal symmetry of the CKM 



matrix V in modulus (i.e., \V as \ ~ \VJ, \V cb \ ~ |F ts | and |FJ ~ \VJ 



In the basis where the flavor eigenstates of three charged leptons are identified 
with their mass eigenstates (i.e., M l — Mj), the Majorana neutrino mass matrix of 
the form 



predicts the /i-r permutation symmetry of the MNSP matrix U with 9 13 = 0° and 
9 23 = 45°; while the mass matrix of the form 

/ a b -b*\ 

M v =\ b c d (7) 
\-b* d c* J 

leads us to the /i-t symmetry of U with 5 = ±90° and 2 3 = 45°. In either of 
the above textures of M„, its entries have certain kinds of linear correlations or 
equalities and thus can be generated from some underlying flavor symmetries. In 
view of the experimental evidence for 9 13 ^ 0° the pattern of M v in Eq. (6) has 
to be modified. For a similar reason, the more reliable and accurate experimental 
knowledge on 9 23 and 5 will be useful for us to identify the effect of /i-r symmetry 
breaking and build more realistic models for lepton mass generation, flavor mixing 
and CP violation. 

3.3. CP and T Violation 

If neutrinos are the Majorana particles, the 3x3 MNSP matrix U contains three 
CP-violating phases 5, p and a. Among them, S determines the strength of CP 
and T violation in neutrino oscillations, because both P{v a — > Vp) — P(p a ^p) 
and P(y a — > Vp) — P(vn — >• v a ) are proportional to the leptonic Jarlskog invariant 
J t = sin# 12 cos6* 12 sin(9 23 cos6* 23 sin# 13 cos 2 9 13 sin 5 in vacuum^]. The phases p 
and a, which have nothing to do with neutrino oscillations, are associated with the 
decay. Note that 8 itself is also of the Majorana nature, although it is usually 
referred to as the Dirac phase: one reason is that 6 may appear in other lepton- 
number-violating processes, even if it can sometimes be arranged not to appear in 
the decay; and the other reason is that <5, p and a are actually entangled with 

one another in the RGE running from one energy scale to another. 

The fact that 9 13 is not strongly suppressed is certainly a good news to the 
experimental attempts towards a final measurement of CP violation in the lepton 
sector. The reason is simply that the strength of leptonic CP violation (i.e., J t ) is 
proportional to sin0 13 . In the qua rk sector one has determined the corresponding 
Jarlskog invariant J ~ 3 X 10~ 5 ^31 anc j attributed its smallness to the strongly 
suppressed values of quark flavor mixing angles (i.e., i? c = ?9 12 ~ 13°, $ 13 ~ 0.2° 
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and $23 — 2.4°). In the lepton sector both 9 12 and 9 23 are large, and thus it is 
possible to achieve a relatively large value of J { if the CP-violating phase S is not 
small either. Taking 9 12 — 34°, 9 13 ~ 9° and 9 23 — 45° as a realistic example of U, 
we arrive at J l ~ 0.036 sin 6, implying that the magnitude of leptonic CP violation 
can actually reach the percent level in neutrino oscillations if S is not strongly 
suppressed. Whether CP violation is significant or not turns out to be an important 
question in lepton physics, especially in neutrino phenomenology. 



3.4. Comparison between the MNSP and CKM Matrices 

The relative sizes of the nine elements of the MNSP matrix U cannot be completely 
fixed unless we have known 9 23 > 45° or 9 23 < 45° as well as the range of 6. With 
the help of the available experimental data and the unitarity of U, we find 

\u el \ > |^ 3 I ~ \u t3 \ > K 2 \ ~ \u T2 \ > \u e2 \ > \u Ml \ ~ \U Tl \ > \U e3 \ , (8) 

where "~" implies that the relative magnitudes of |J7 Mi | and \U Ti \ (for i = 1,2,3) 

remain undetermined at present. In comparison, t he nin e elements of the CKM 

170171 1 

matrix V are known to have the following hierarchy 1 1 \ 

\V tb \ > \V ud \ > \VJ » \V US \ > \VJ » \V cb \ > \V ts \ » \V td \ > \V ub \ . (9) 

There is a striking similarity between the quark and lepton flavor mixing matrices: 
the smallest elements of both V and U appear in their respective top-right corners. 

In the history of flavor physics it took quite a long time to measure the four inde- 
pendent parameters of V, but the experimental development had a clear roadmap: 

*12 (o r Ksl) => ^23 (or \V cb \) 13 (or \V ub \) 5 (quark) . (10) 

Namely, the observation of the largest mixing angle d 12 was the first step, the 
determination of the smallest mixing angle i? 13 was an important turning point, and 
then the quark flavor physics entered an era of precision measurements in which CP 
violation could be explored and new physics could be searched for. Interestingly and 
hopefully, the lepton flavor physics is repeating the same story: 

#23 (or l^sl) 012 (or \U e2 \) =^> 9 13 (or \U e3 \) S (lepton) , (11) 

where 9 23 is the largest and 9 13 is the smallest. The observation of 9 13 in the Daya 
Bay experiment is paving the way for future experiments to study leptonic CP 
violation and to look for possible new physics (e.g., whether the 3x3 MNSP matrix 
U is exactly unitary or not LL ^), in particular through the measurements of neutrino 
oscillations for different sources of neutrino beams. The Majorana nature of three 
massive neutrinos and their other two CP-violating phases (i.e., p and a) can also 
be probed in the new era of neutrino physics. 
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4. POSSIBLE LEPTON FLAVOR STRUCTURES 
4.1. Two Phenomenological Strategies 

The MNSP matrix U actually describes a fundamental mismatch between the flavor 
and mass eigenstates of six leptons, or a mismatch between diagonalizations of the 
charged-lepton mass matrix M l and the effective neutrino mass matrix M v in a given 
model, no matter whet her the origin of neutrino masses is attributed to the seesaw 
mechanisms or not ^1 Assuming massive neutrinos to be the Majorana particles, 
we may simply write out the leptonic mass terms as 




Anass = (e' /*' r'V M, U' + - (y e u„ v T ) M v i£ + h.c. , (12) 



where "/" stands for the flavor eigenstates of charged leptons, "c" denotes the charge- 
conjugated neutrino fields, and M v is symmetric. By using the unitary matrices ; , 
0[ and O v , one can diagonalize M l and M v through the transformations 0\M i O\ = 
M l = Diag{m e , m„, m T } and 0\M v O* = M v = Diag{m 1 , m 2 , rn 3 }, respectively. 
Then one arrives at the lepton mass terms in terms of the mass eigenstates: 

1 




C = (e A* t) l Mi | n ) + - ( Vl v 2 v 3 ) L M u \v%\ + h.c. . (13) 

a V^Iy R 

Extending this basis transformation to the standard charged-current interactions, 
we immediately obtain 

-£ cc = -^ (e M r) L 7^[^] ^" + h.c. , (14) 

in which U — 0\O v . The above treatment is most general at a given energy scale 
(e.g., the electroweak scale). There are two different strategies of phenomenologically 
understanding the structure of the MNSP matrix ^1 

(1) The mixing angles of U are associated with the lepton mass ratios. The 
structure of lepton flavor mixing is directly determined by the structures of O l 
and O v . Since these two unitary matrices are used to diagonalize M t and M v , 
respectively, their structures are governed by those of M l and M v , whose eigenvalues 
are the physical lepton masses. Therefore, we anticipate that the dimensionless flavor 
mixing angles of U should be certain kinds of functions whose variables include four 
independent mass ratios of three charged leptons and three neutrinos. Namely, 

•J , (15) 

where the Greek subscripts denote the charged leptons, the Latin subscripts stand 
for the neutrinos, and "• • • " implies other dimensionless parameters originating 
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from the lepton mass matrices. Such an expectation has proved valid in the quark 
sector to explain why the relation sin"i? c ~ \J m d/ m s works quite well and how 
the hierarchical structure of the CKM matrix V is related to the st rong hierarchies 
of quark masses (i.e., m u <C m c <C m t and m d <C m s <C m b ) ^1 As for the 
phenomenon of lepton flavor mixing, it is apparently difficult to link two large 
mixing angles 9 12 and 9 23 to m e /m„ ~ 4.7 x 10 -3 and m„/m T ~ 5.9 x 10 -2 IZSIlS 
Hence one may consider to ascribe the largeness of 9 12 and 9 23 to a weak hierarchy 
of three neutrino masses, such as the conjecture tan# 12 ~ yjm 1 /m 2 1 78 | 79 | 

To establish a direct relation between 9^ and lepton mass ratios, one has to 
specify the textures of M l and M v by allowing some of their elements to va nish or 
to be vanishingly small. A typical example of this kind is the Fritzsch ansatz I 80 | 81 | ; 

/0 x 0\ 

M l>v = x x , (16) 
\0 x x) 

which is able to account for the present neutrino oscillation data to an acceptable 
degree of accuracy (e.g., sin# 23 ~ J m^/m T + ^Jm 2 jm 3 ~ 0.65) l&2J 83 | 84 | Another 
well-known an d viable exa mple is the two-zero textures of M v in the basis where 
M l is diagonal ESESIBl lSBl Note that the texture zeros of a fermion mass matrix dy- 
namically mean that the corresponding matrix elements are sufficiently suppressed 
as compared with their neighboring counterparts, and they can be derived from a 
certain flavor symmetry in a giv en theoretical framework (e.g. , w ith the help of the 
Froggatt-Nielson mechanism ISHI or discrete flavor symmetries I^QI) . 

(2) The lepton flavor mixing matrix U consists of a constant leading term U 
and a small perturbation term AU. In f act, U h as been conjectured to have the 
following structure for a quite long time | 73 | 91 | 92 [ 

U = (U + AU) P u , (17) 

in which the leading term U is a constant matrix responsible for two larger mixing 
angles 9 12 and 9 23 , and the next-to- leading term AU is a perturbation responsible 
for both the smallest mixing angle 9 13 and the Dirac CP-violating phase 5. So far 
a lot of flavor symmetries have been brought into exercise to derive U , while AU 
might originate from either an explicit flavor symmetry breaking scenario or some 
finite quantum corrections at a given energy scale or from a superhigh-energy scale 
to the electroweak scale. 

In this case the MNSP matrix U is approximately a constant matrix whose mix- 
ing angles are independent of the lepton mass ratios. This conjecture is actually in 
conflict with the conjecture made in the first strategy. The reason for this "conflict" 
is rather simple: the assumed structures of lepton flavor mixing in Eqs. (15) and 
(17) correspond to two different structures of lepton mass matrices. As we have 
pointed out above, the direct dependence of 9^ on m a jm^ and m k /m l is usually a 
direct consequence of the texture zeros of M ; and (or) M v . In contrast, a constant 
flavor mixing pattern U may arise from some special textures of M l and (or) M v 
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whose entri es have certain kind s of linear correlations or equalities. For instance, 



assures O v to be of the tri-bimaximal mixing pattern to be discussed in section 
4.2. This neutrino mass matrix has no zero entries, but its nine elements satisfy 



the sum rules {M v ) a + {M v ) i2 + (M u ) i3 = and (MJ y + {M„) 2j + (M„) Sj = 



(for i,j = 1, 2, 3). Such correlative relations are similar to those texture zeros in the 
sense that both of them may reduce the number of free parameters associated with 
lcpton mass matrices, making some predictions for the lepton flavor mixing angles 
technically possible. 

In short, one may try to understand the structure of the MNSP matrix U by 
following two phenomenological strategies: 

(1) to explore possible relations between the flavor mixing angles and the lepton 
mass ratios; 

(2) to investigate possible constant patterns of lepton flavor mixing as the leading- 
order effects. 

The first possibility points to some vanishing (or vanishingly small) entries of M ; 
and M u , while the second possibility indicates some equalities or linear correlations 
among the entries of M l or M v . In both cases the underlying flavor symmetries 
play a crucial role in deriving the structures of lepton mass matrices which finally 
determine the structure of lepton flavor mixing. Of course, how to pin down the 
correct flavor symmetries remains an open question. 

4.2. Five Typical Patterns of U 

It is well known that the special textures of M l and M„ like that in Eq. (18 ) can 
easily be derived from certain discrete flavor symmetries (e.g., A 4 or 5 l 4 )| 93|94 l That 
is why Eq. (17) formally summarizes a large class of lepton flavor mixing patterns 
in which the leading terms are constant matrices originating from some underlying 
flavor symmetries. The fact that 13 is not very small poses a meaningful question 
to us today: can this mixing angle naturally be generated from the perturbation 
matrix At/? The answer to this question is certainly dependent upon the form of 
U in the flavor symmetry limit. Here we reexamine five typical patterns of U in 
order to get a feeling of the respective structures of AU which can be constrained 
by current experimental data on neutrino oscillations. 

For the sake of simplicity, we typically take 9 12 ~ 34°, 9 13 ~ 9° and 9 23 — 45° 
as our inputs to fix the primary structure of the MNSP matrix U. Then we have 



the texture 



161162 163 64 65 l66l67l 




(18) 




(19) 
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It makes sense to compare a constant mixing pattern U with the observed pattern 
of U in Eq. (19), such that one may estimate the structure of the corresponding 
perturbation matrix AU. Let us consider five well-known patterns of U in the 
following for illustration. ^^^^^^ 
(1) The democratic mixing pattern of lepton flavors [Z 3 | 91 | 92 |. 



Co=|-£ 3* 5§h (20) 

3(0) _ ko d°) - no n „A o(°) - arctan( v / 2) 54.7° 
in the standard parametrization as given in Eq. (2). With the help of Eq. (19), we 
immediately obtain the form of AU = UP} — U as follows: 

0.112 -0.155 0.156e-^ 

AU = | 0.013- 0.092e l5 0.178 - 0.062e l5 -0.118 | . (21 i 

-0.182 - 0.092e 4<5 -0.009 - 0.062e l5 0.121 

One can see that the magnitude of each matrix element of AU is of 0(0.1), implying 
that the realistic pattern of U might result from a democratic perturbation to U Q 
(i.e., the nine entries of AU are all proportional to a c ommon small parameter). 
(2) The bimaximal mixing pattern of lepton flavors E 5 | 96 [ 



Uo = -3 5 3s > ( 22 ) 




which has 9^ = 45°, ofj = 0° and 9^ = 45° in the standard parametrization. 
Comparing Eq. (22) with Eq. (19), we obtain the perturbation matrix 

0.112 -0.155 0.156e~ J,5N 

AU = | 0.105- 0.092e l5 0.086 - 0.062e l5 -0.009 | . (23) 
-0.105 - 0.092e l5 -0.086 - 0.062e l5 -0.009 

We see that the matrix elements (AU) „ 3 and (AU) t3 are highly suppressed. In 
other words, the initially maximal angle 9^ receives the minimal correction, which 
is much smaller than the one received by the initially minimal angle 9^ . Such a 
situation is more or less unnatural, at least from a point o f view of mod el building. 
(3) The tri-bimaximal mixing pattern of lepton flavors IS 7 | 98 | 99 | l IIIl 



U = 



I V2 J_ n 

' V3 \/3 U 

1 1 1 

V6 %/3 \/2 

\ J_ 1 J_ 



(24) 



whose three mixing angles are 9^ = arctan(l/v / 2) — 35.3°, 9^ = 0° and 9^ — 
45° in the standard parametrization. In a similar way we get the corresponding 
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perturbation matrix 



0.003 

AU = 1 0.013- 0.092e 115 
-0.013 - 0.092e i5 



-0.025 0.156e- l5N 
0.009 - 0.062e l5 -0.009 
-0.009 - 0.062e il5 -0.009 



(25) 



It is quite obvious that (AU) el , (AU) e2 , (AU) 3 and (AU) t3 are highly suppressed. 



So two initially large angles 9^ and 0^' are only slightly modified by the pertur- 
bation effects, but the initially minimal angle 9^ receives the maximal correction. 



3 (0) 



(4) The golden-ratio mixing pattern of Icpton flavors 



110111021 



V2 



V2 



a/5-V5 
1 



1 



1 1 



1 

v/2 
1 



(26) 



which has 9 



(0) 
12 



arctan[2/(l + a/5)] ^ 31.7°, 0$ = 



and 4? = 45° in the 



standard parametrization. In this case the perturbation matrix AU turns out to be 

0.026 0.156e~ l,5N 
-0.016 - 0.062e i<5 -0.009 1 . (27) 



AU = 



-0.032 
-0.023 - 0.092e 4<5 
0.023- 0.092e l5 0.016 



0.062e 



-0.009 



Similar to the tri-bimaximal mixing pattern, two initially large angles of the golden- 
ratio mixing pattern are only slightly corrected, but the initially minimal angle 9^ 
is significantly modified by the same perturbation. ^^^^^^ 
(5) The hexagonal mixing pattern of lepton flavors H 03 1 1 04 | 1 5l 



(28) 




(o) _ 



30° 



9\l ] = 0° and 6» 2 3 ; = 45° in the standard 



3(0) _ 



whose mixing angles are 9\ 
parametrization. In this case we obtain the perturbation matrix 

-0.047 0.052 0.156e- 4,5N 

Af = I -0.041 - Q.Q92e lS -0.026 - 0.062e l5 -0.009 
0.041 - 0.092e l5 0.026 - 0.062e i<5 -0.009 



(29) 



This result is quite analogous to the one obtained in Eq. (25) or Eq. (27), simply 
because the patterns of U in these three cases are quite similar. 

Now let us summarize some useful lessons that we can directly learn from the 
above five typical examples of U. 

• To accommodate the observed value of 9 13 in a generic flavor mixing structure 
U = (U + AU) P v , one has to choose a proper constant mixing pattern U Q and 
adjust its perturbation matrix AU. The phenomenological criterion to do so 
is two-fold: on the one hand, U Q should easily be derived from a certain flavor 
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symmetry; on the other hand, AU should have a natural structure which can 
easily be accounted for by either the flavor symmetry breaking or quantum 
corrections (or both of them). 

The common feature of the above five patterns of U is apparently (U ) e3 = 
(or equivalently, 8^ = 0°), implying that a relatively large perturbation is 
required for generating 13 ~ 9°. In this case, the closer 0$ and #23 arc to the 
observed values of 9 12 and 9 23 , the more unnatural the structure of AU seems 
to be. The tri-bimaximal mixing pattern given in Eq. (24), which is currently 
the most popular ansatz for model building b ased on certain flavor symmetries, 
suffers from this unnaturalness in particular EE3 i n this sense we argue that 
the democratic mixing pattern in Eq. (29) might be more natural and deserve 
some more attention. 

One may certainly consider some possible patterns of U a which can predict a 
finite value of 6^ in the vicinity of the experimental value of 9 13 . In this case 
the three mixing angles of U may receive comparably small corrections from 
the perturbation matrix AU, and thus the naturalness criterion can be satisfied. 
For example, the following two patterns of U belong to this category and have 
been discussed in th e literature EnZHDSI. One of them is the so-called correlative 
mixing pattern tlQQl 



u =\-js- -W" 5 4g - -W* ^ (30) 




with c, = cos#„, = (V2 + l)/VG and s„ = sm0* = (V2— l)/V&, which predicts 
9$ = arctan(l/V2) -35.3°, 9 { 23 } = 45° and 0$ = 9 { 2 ° 3 } - 9 { $ ~ 9.7°. The three 
mixing angles in this constant scenario satisfy two interesting sum rules: 

#12 + #13 = ^23 > 

^12+ #13 + ^3 = 90°, (31) 

which are geometrically illust rated in Fig. 3. The other pattern of U Q is the 
tetra-maximal mixing pattern HHHI 



u 



( 2+V2 1 2-\/2 \ 

4 2 4 

4 4 2 44" 1 " 4 

. _ ^2 _ <V2-l) 1 iV2 VI i(V2+l) . 

\4 4 2'44 4/ 



(32) 



which predicts 9 { ° 2 ] = arctan(2 - \/2) ~ 30.4°, 9 { 23 } = 45° and 9^ = arcsin[(2 - 
V2)/4] ~ 8.4°. Of course, whether such constant mixing patterns can easily be 
derived from some underlying flavor symmetries remains an open question. 

In short, today's model building has to take the challenge caused by the reasonably 
large value of 9 13 . 
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6 V1 ~ 35.3° 




Fig. 3. A geometrical description of the sum rules 9 12 + 13 = 6 2 3 anc ' ^12 + ^13 + ^23 = 90° for 
the correlative neutrino mixing pattern in terms of the inner angles of the right triangle AABC. 

Note that the RGE running effects or finite quantum corrections are not easy 
to generate 9 13 ~ 9° from 8^ = 0°. unl ess the seesaw threshold effects or other 
extreme conditions are taken into accounti 110 ! 111 ! 112 ! 113 ! 114 ! 115 ! 116 ! 117 ! 118 ! 119 ! 
One may therefore consi der a pattern of U with nonzero 9^ , such as the tetra- 
maximal mixing pattern EQ] or ^- ne correlative mixing pattern H2I] as a starting 
point of view to calculate the radiative corrections before confronting it with current 
experimental data. We shall elaborate on this point in detail in section 5. 

4.3. The Minimal Perturbation to U 

Note that the perturbation matrix AU in Eq. (17) is in general a sum of all possible 
perturbations to the constant flavor mixing matrix U . From the point of view of 
model building, it is helpful to single out a viable AU whose form is as simple as 
possible. To do so, let us reexpress Eq. (17) in the following manner: 

U = (U + AU) P v = U Q (1 + AU') P v = (l + AUl) U (1 + AC74) P v , (33) 

where AU = U Q AU' = AU{U Q + U AU{ t + AU[U Q AU^ holds, and it satisfies the 
condition U AW + AUUl + AUAW = as a result of the unitarity of U itself. 
Therefore, one may achieve a viable but minimal perturbation to U by switching off 
A/7l (or AU^) and adjusting AU^ (or AU{J to its simplest form which is allowed 
by current experimental data. Such a treatment is actually equivalent to multiplying 
U by a unitary perturbation matrix, which may more or less deviate from the unit 
matrix 1, from either its l eft-hand side or its right-hand side. The first example of 
this kind was given before H 3 | 91 | 92 | f Qr ^ ne democratic mixing pattern, and its AU 
was mainly responsible for the generation of nonzero 9 13 and 6. 

Here we concentrate on the typical patterns of U discussed above and outline 
the main ideas of choosing the minimal perturbations to them. 

• If U predicts 8 { 2 ° 3 } = 45° and 9 { ° 3 ] = 0° together with 6$ > 34° (the best-fit 
value based on current neutrino oscillation data , then the simplest way to 
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generate a relatively large 9 13 , keep 9 23 = 9 23 = 45° unchanged and correct 
9^2 to a slightly smaller value is to choose a complex (2, 3) rotation matrix as 
the perturbation matrix: 



'10 
AU' = I cos6> isinf 
\0 i sin 9 cos 9 



AU' 



-±sin 2 l 
v0 ism 9 



1 sm f 
-5 sin 2 1 



(34) 



where 9 is a small angle to trigger the perturbation effect. The most striking 
example in this category is to take t/n t o be the tri-bimaximal mixing pattern 
given in Eq. (24). The result isHHUSl 

y/2 1 



V3 V3^° W 

r = '-To 73 COS0 + 7f 
7e -75 COS0 + 71 



smt 



smf 



sine 



1 



cos ( 



cos ( 



^sint? | P„ 
73 SiE 



(35) 



which predicts 
1 



(l - 2tan 2 t9 13 ) , sin 2 



1 



sin 2 9 , 9 23 = 45° , 6 = 90° (36) 



in the standard parametrization. Note that the obtained correlation between 
9 12 and 9 13 is especially interesting because it leads us to 9 12 — > 34° when 
9°, consistent with the present experimental data. II 9 23 is allowed to 



'13 



slightly deviate from 9 2 ° 3 = 45°, then one may simply make the replacement 



i ->■ e ld in Eq. (35). 

• If U predicts 9% = 45° and 6$ = 0° together with 6$ < 34°, then the 
most economical way to generate a relatively large 9 13 , keep 9 23 — 9^ 23 = 45° 
unchanged and correct 9f 2 to a slightly larger value is to choose a complex 
(1,3) rotation matrix as the perturbation matrix: 



AU' 



' cos 9 i sin 0\ 

1 
K isin9 cos( 



or AU' 



■ - \ sin 2 i 


i sin 9 



% sm f 




(37) 



-±sin 2 i 



Taking U to be the golden-ratio mixing pattern in Eq. (26), we immediately 
arrive at 

V2 is/% 



U = 



_V2_ 

1 cos 9 



COS( 



1 



y/2' 



coat 



whose predictions include 9 23 
2 



sm f — 



\/5+\/5 

1 1 

i i 



\/5— v/5 



sm f 



45°, (5 = 90°, and 



COS ( 

cos( 



, sm (/ 1 

V5+V5 / 



P„ , (38) 



sm 



'12 



1 +tan z 



sm 



(39) 



5 + 75 v ' id 5-^5 

in the standard parametrization of U. In this case the correlation between 9 12 
and 9 13 leads to 9 12 — > 32° when 9 13 — > 9°, compatible with the experimental 
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data. Again, the replacement i — > e lS in Eq. (38) allows one to obtain a some- 
what more flexible value of #23 which may slightly deviate from ^3 = 45°. 
• If U Q is quite far away from the realistic MNSP matrix U, one has to consider 
a somewhat complicated perturbation matrix including two rotation angles. In 
the neglect of CP violation, for instance, we may consider 

/ c' 12 s'u \ 
1 + AU' = s[ 2 ^ 3 c[ 2 c' 23 s' 23 , (40) 

\ s 12 s 23 c 12 s 23 — c 23/ 

where = cos 9'^ and s'^ = sin 6?^- (for ij = 12,23). However, we hope that 
the resulting structure of U still allows us to obtain one or two predictions, in 
particular for the mixing angle 9 13 . A simple example of this kind has been given 
before 1 1 by taking U to be the democratic mixing pattern, and it predicts an 
interesting relationship between 9 13 and 9 23 m the standard parametrization: 

, \/2-tan0 2 3 . . 

am0 13 = . (41) 

J 5- 2\/2 tan 9 23 + 4 tan 2 9 23 

Typically taking 9 2 $ — 45°, we can arrive at 9 13 ~ 9.6° \i \ s certainly easy 
to accommodate a CP-violating phase in AC/', although its form might not be 
really minimal anymore. 

For those constant flav or m ixing patterns with 9v^ ^ 0° from the very beginning, 



such as the correlative and tetra-maximal^®! mixing scenarios given in Eqs. 
(30) and (32), the similar minimal perturbations can be introduced in order to 
make the resulting MNSP matrix U fit the experimental data to a better degree of 
accuracy. 

It should be noted that the above discussions about possible patterns of AU 
(or AC/') with respect to those of U are purely phenomenological. From the point 
of view of model building, it is more meaningful to consider the textures of lepton 
mass matrices 

M l = Af/ 0) + AM; , M v = + AM V , (42) 

where M ; (0) and M^ 0) can be obtained in the limit of certain flavor symmetries, and 
their special structures allow us to achieve a constant flavor mixing pattern U . The 
perturbation matrices AM l and AM V play an important role in transforming C7 
into the realistic MNSP matrix U, and thus their textures should be determined 
in a simple way and with a good reason. The connection between AM l v and AU 
(or AC/') depends on the details of a lepton flavor model and may not be very 
transparent in most cases. In the basis where M l is real and positive, however, 
AM„ can be formally expressed as 

AM V = (U + AU) ~M V (U + AU) T - U ~Mf C/ T , (43) 
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where ~M V = P V M V P^ and M { ° ] = P' V M' V P'J , M' v = Diag{mi, m 2» m 3> and K = 
Diag{e lp ', e l<J , 1}. Here m- (for i = 1,2,3) denote the eigenvalues of M„ in the 
symmetry limit, while p' and a 1 stand for the Majorana phases in the same limit. 
It is therefore possible, at least in principle, to fix the structure of AM V with the 
help of a certain flavor symmetry and current experimental data. 



5. UNSUPPRESSED 6» 13 AND RGE RUNNING EFFECTS 

The RGE running effects of the neutrino flavor parameters have been discussed 

in many pap ers 

lll()lllllll2lllHlll4lll5lll6lll7lll8lliail25ll26ll27ll28ll29ll3()ll31llM2ll^Hll^4l It 

is known that large radiative corrections to those parameters are possible, especially 
when the neutrino masses are nearly degenerate or the value of tan /3 is sufficiently 
big in the MSSM. The fact that 6 13 is not as small as previously expected motivates 
one to reconsider how it can be generated at the tree level or by quantum correc- 
tions. Some studies in this regard have recently been done to look at the impacts of 
a relatively large 9 13 on the running behaviors of the other two mixing angles and 
the CP-violating phases tt^H. 



5.1. Approximate RGEs 

The masses of the Majorana neutrinos are believed to be attributed to some un- 
derlying ne w physics at a superhigh-energy scale A (e.g., via the canonical seesaw 
mechanism | 20 | 2 1 | 22 | 23 | 24 ]^ -g^ j, mc j f new phy S i cs can all point to the unique 
dimension-5 Weinberg operator for the neutrino masses in an effective field theory 
after the corresponding heavy degrees of freedom are integrated out HHH. I n the 
MSSM, such a dimension-5 operator reads 

= \kH 2 ■ « • Hjit + h.c. , (44) 



where A denotes the cutoff scale, £ L stands for the left-handed lepton doublet, 
H 2 is one of the MSSM Higgs doublets, and K represents the effective neutrino 
coupling matrix. One may obtain the effective Majorana neutrino mass matrix M v = 
kv 2 tan 2 /3/(l+tan 2 0) after spontaneous gauge symmetry breaking. The cutoff scale 
A implies the scale of new physics, such as the mass scale of the heavy Majorana 
neutrinos in the canonical seesaw mechanism. The evolution of K from A down to 
the electroweak scale A EW is formally independent of any details of the relevent 
model from which n is derived. Below A the scale dependence of n is described by 



167T — = a M n + 
at 



(r i y i t ) K + K (y ; y/) 



(45) 



at the one-loop level in the MSSM H25l l^ | 128 | w h ere a M rj ~l.2gl - 6g 2 + 6yf. 
One may use Eq. (45) to derive the explicit RGEs of the three neutrino masses 



3 tnc explicit riuiis oi tin 
[TT5IT 14 125 126 128TT5T1 



and six flavor mixing parameters l 11 ° l 11 ^ l lz, - > l lz '- ) l lzo l 1J1 [ Given an approximate 
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mass degeneracy of the three neutrinos together with the standard parametrization 
of the MNSP matrix U in Eq. (2), the RGEs of m, (for i = 1,2, 3) turn out to be 

dTTT/^ TTi-^ r 2/22 /r\ / 2 \\1 

~dT ~ 167T 2 \- ayl + Vt v s i2 s 23 - 2 c s c 12 c 23 s 12 s 23 s 13 + C(s 13 )jJ , 
dm 2 m 2 



dt 16tt 2 
dm 3 m 3 



[a M + 2y 2 (c 2 2 s 23 + 2c,5C 12 c 23 s 12 s 23 s 13 + Q(s 2 13 ))] , 



2y 2 c 2 23 + 0(s 2 13 )] 



dt ~ 16tt~ L " m 
The RGEs of 9 tJ (for ij = 12, 23, 13) are found to be 

d0, 



(46) 



y 12 
dt 



( C 12 s 12 s 23 — cos 2^i 2 C 23 S 13 C 5 ) C( p _ CT ) 



66. 



23 



dt 



Am 32 

2 2 
2 
1 



+ c 23 s 13 s <5 s (p-o-) 
i 3 _c 23 s 23 s 13 ( s 12 c (5+p) c p + c 12 c (<5+<r) c cr) + C( s 13) j i 



j/r "1 f / 2 2 . 2 2\ 



-C 12 S 12 C 23 S 13 I C( 5+ 2p) C (<5+2cr) 



0(4) 



d0 13 y 2 to 2 



dt ~ 87r 2 Am| 2 L23Ll3 



C 12 S 12 S 23 I C ((5+2p) C (<5+2<r) 



_ 2c 23 s 13 ^c 2 2 c 2 5+p j + s 2 2 c 2 5+cr ^ + 0(s 2 3 ) 



(47) 



in which c x = cos X and s.,. = sinx (for x = 5, p, u, p—cr, 5+p, 6+cr, S+2p, S+2<r). 
The RGEs of the three CP- violating phases 5, p and a can be written as 



dt ~ 4?r 2 {Am^ 



9 "23 



cos 26* 12 c 23 s 13 c 5 



'12 a 12 



C 23 S 13 S <5 



12 s 12 



>(p-<0 



2 C 12 S 12 C 23 S 23 ( S (t5+2p) S (c5+ 2 - 
+ ( c (5-p) s (5-p) s 12 + c («5- -) s ((5- ( t) c 12) cos 2^ 23 S 13 

+ (cpS p c 2 2 + c a s a s\ 2 ) c 2 23 s 13 + 0(s 2 3 )] } , 



dp _ y 2 T f m| 2 
dt ~ 4tt 2 \Am^ 23Cl2 



S 23 — 



cos 2# 12 c 23 s 13 c (5 



'(/>-*) 



c 23 s 13 s <5 



C-i nS 



12 S 12 



»(p-<0 
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( c (<5-p) s (<5-p) s 12 + c (8-cr) s (S-a) c 12^ cos2 # 2 3 s lS 

cos20 12 c 2 3S 13 c 5 



0{s\ 3 ) 



2 S 23 S 12 



C-inS 



12°12 



-{p-a) 



+ - 



c 23 s 13 S i5 
c 12 s 12 



• (p _ ff) +o( s ? 3 ) 



s (5-<t) c i 2 ) cos26» 23 s 13 + C(s 



13/ 



In addition, the RGE of J' is obtained as follows: 



dt 



J 



!Jr 



8tt 2 A 



~1T~ C0S 2^12 S 23 — cos2 2^12 C 23 S 23 C 13 S 13 C <5 S I 5 



'21 



A, 



13, 



'32 



(48) 



(49) 



The running behaviors of the three mixing angles and three CP-violating phases for 
a very small 9 13 can be very different from those for a relative large 9 13 . In particular, 
the CP-violating phases play a crucial role in the RGEs. Let us elaborate on this 
point in the following. 



5.2. Flavor Mixing Angles and CP-violating Phases 

(1) The running behaviors of the three mixing angles. Eq. (47) shows that the RGE 
running behaviors of the three neutrino mixing angles are strongly dependent on the 
three CP-violating phases. As for the Majorana neutrinos, the radiative corrections 
to the three mixing angles can be adjusted by choosing different values of the CP- 
violating phases 8, p and a. A numerical analysis has been carried out to look 
at their numerical evolution to A FS via the RGEs in the MSSM with tan/3 = 10 
(denoted as "MSSM10" for short) or tan/3 = 50 (denoted as "MSSM50" for short), 
in which 8 12 = 34°, 9 23 = 46°, # 13 = 9° are taken as the typical input and the 
three CP-violating phases are freely adjusted at A EW G2H. The main results are 
summarized in Table 2, where the upper (lower) lines show the possible ranges of 
three mixing angles at A FS for the normal (inverted) neutrino mass hierarchy. 

2 2 

The RGE of 9 12 is dominated by the term — — ^-c^s^s^cf implying 

8n z Ato 2 i 

that the magnitude of the radiative correction to 9 12 depends strongly on the phase 
difference (p — a). Hence 9 12 is most sensitive to the RGE effect when p ~ a holds. 
Running from A FS ~ 10 14 GeV down to A EW , the mixing angles 9 23 and 9 13 receive 
less significant radiative corrections in the MSSM10 case, as shown in Table 2. While 
in the MSSM50 case 9 23 and 9 13 may also receive significant radiative corrections 
if three CP-violating phases are well turned. 
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Table 2. The radiative corrections to the three mixing angles from 
A EW ~ 10 2 GeV to A F g ~ 10 14 GeV in the MSSM with tan ,8 = 10 
or 50, where the three CP-violating phases are freely adjusted. 



Parameter 


Input (A EW ) 


Output (A FS ) 










MSSM10 


MSSM50 






34.0° 


7° to 55° 
2° to 31° 


0.5° to 62° 
2° to 45° 


(NH) 
(IH) 


9 23 


46.0° 


47° to 48.5° 
43.5° to 45° 


7.5° to 45.5° 
46.5° to 89° 


(NH) 
(IH) 




9.0° 


8.5° to 10° 
7.5° to 9.5° 


2° to 23° 
7.5° to 82° 


(NH) 
(IH) 



We have seen that the values of the three CP- violating phases are crucial for the 
evolution of the three mixing angles. A very special case is (p — a) ~ ±90°, which 
leads us to 
d0 12 y 2 T m\ 



( 2 i 2 

n 2 " C 23 S 23 S 13 \ S 12 c (S+p) C p + C 12 s (8+p) S p 
I ^ 

d6 *23 _ Vr m l 



dt 4ir 2 Am 32 



2 ™2 

(2 2,2 2\ 
23 \ s \2 c p + c 12 s p) c 12 s 12 c 23 S 13 C (<5+2p) 



dt ~ 47r 2 Ami 2 C23 
d0 13 yi ^ 



" c 23 c 13 
'32 



12 s 23 C (5+2p) C 23 S 13 ( C 12 c («5+p) + S 12 s (<5+p)) -(50) 



dt 

Note that the term proportional to mf/Am 2 ,! in the RGE of 12 in Eq. (47) is 
suppressed by cos(p — a) ~ in this special case, and thus it has been omitted 
from Eq. (50). The three mixing angles may therefore receive comparably small 
radiative corrections for a modest value of tan/3 (e.g., in the MSSM10 case). This 
observation was not noticed in th e literature simply because 9 13 used to be assumed 
to be very small | 125 | 126 | 127 | 131 | if tan p [ s sufficiently large (e.g., in the MSSM50 
case), however, the phase difference (p — a) will be able to quickly run away from 
its initial value (p— cr) ~ ±90° due to the significant radiative corrections, implying 
that Eq. (50) is no more a good approximation of Eq. (47). 

Here let us consider a typical example of this special case — the correlative 
neutrino mixing pattern with 9 12 ~ 35.3°, 9 23 = 45° and 9 13 ~ 9.7°Q2H3 Compared 
with the best- fit values of the three mixing angles at A EW , the three mixing angles 
in this correlative mixing pattern at A FS have to receive comparably small radiative 
corrections during their RGE evolution. As we have mentioned in the last paragraph, 
this requirement can easily be achieved in the MSSM10 case provided the condition 
(p — a) ~ ±90° is satisfied for a nearly degenerate neutrino mass spectrum. Such a 
condition is unnecessary if the neutrino mass spectrum has a strong hierarchy. The 
numerical results are presented in Table 3, where 5 = —68°, p = 13° and a = 115° 
are input at A FS - Then 12 = 34.52°, 9 23 = 45.98° and 13 = 8.83° are obtained at 
Aew after the radiative corrections. 

(2) The radiative generation of the CP-violating phase s. It is well known that 
one CP-violating phase can be generated from another simply because they are 
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Table 3. Radiative corrections to the correlative neutrino mixing 



pattern from A pg ~ 10 14 
Parameter 


GeV to A EW ~ 10 2 
Input (A FS ) 


GeV in the MSSM10. 
Output (A EW ) 


m x (eV) 


0.227 


0.200 


Am| x (10" 5 eV 2 ) 


15.72 


7.59 


Am 2 ^ (10" 3 eV 2 ) 


3.19 


2.40 


e \i 


35.3° 


34.52° 


e 23 


45° 


45.98° 


9 13 


9.7° 


8.83° 


S 


-68° 


-80.88° 


P 


13° 


19.64° 


a 


115° 


118.03° 


8 + p + a 


60° 


56.79° 



entangled in the RGEs. An especially interesting example is the Dirac phase 5, which 
measures the strength of CP violation in neutrino oscillations at the electroweak 
scale, can be radiatively generated from the nonzero Majorana phases p and a at 
a supcrhigh-energy scale. If 9 13 is very small, however, the running of 6 can be 
significantly enhanced by the terms that are inversely proportional to sin 9 13 . In the 
MSSM10 case it has b een found that even S — 90° can be radiatively generated if 
9 13 ~ 1° is taken USSl. Given 9 13 ~ 9° at A EW , it is found that -30° < S < 30° 
at A E w may result from 6 = 0° at A FS inthe MSSM10 case. In the MSSM50 case 
even \8\ ~ 90° can be obtained at Aew l^fl 

(3) The running of the sum 5 + p + a. Eq. (48) leads us to the RGE of the sum 
of the three CP-violating phases: 

d rx j_ „ ^ ~ m i 1 

— {5 + p + a) - 



AV 47r 2 Am| 2 s 13 



'13 



2 L 12°'12 L 23''23 y*(5+2p) "(S+2a) / 
~ ( c (S-p) s (8-p) s 12 + c {S-a) S {5-a) C l2) COs29 23Sl 

+ {c p s p cl 2 +c aSa sl 2 )c 2 23Sl3 + 0{s 2 13 )] . (51) 

Since the value of 9 13 is not small, the RGE running effect on (8 + p + a) is expected 
to be insignificant. In other words, the sum of the three CP- violating phases may 
approximately keep unchanged during the RGE evolution in the SM or MSSM with 
a modest tan j3. The numerical analysis shows that [6 + p + a) changes less than 
4° when running from A FS down to A EW in the MSSM10 case. The stability of 
(5 + p + a) aga inst the radiative corrections is quite impressive, unless tan/3 is 
sufficiently large 11211 

To summarize, in order to obtain a phenomenologically-favored neutrino mixing 
pattern at the electroweak scale A EW , the radiative corrections should be carefully 
examined for those mixing patterns at a superhigh-energy scale A FS which might 
result from a certain flavor symmetry. In the MSSM10 case the values of 6*23 and 9 13 
predicted at A FS are always close to their running values at A EW , while the value 
of 9 12 at A FS can be somewhat smaller or larger than its running value at A EW . In 
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the MSSM50 case the allowed ranges of the three mixing angles at A FS can be quite 
wide, as we have discussed above. However, a crucial point is that a given flavor 
symmetry model should be able to predict the appropriate CP-violating phases 
at A FS in order to obtain the appropriate mixing angles at A EW after the RGE 
evolution. In general, it is possible to generate 6> 13 ~ 9° at A EW from 8 13 ~ 0° at A Fg 
through the radiative corrections, in particular when some new d egrees of free dom 
or nontrivial running effects (such as the seesaw threshold effects ^ 13 | 114 |131|^ are 
taken into account or the three CP-violating phases are fine-turned. Therefore, we 
argue that it seems more natural for a specific flavor symmetry model to predict a 
relatively large 9 13 at Afs- 

A measurement of the Dirac phase 5 in the forthcoming long-baseline neutrino 
oscillation experiments and any experimental information about the Majorana CP- 
violating phases p and a are extremely important, so as to distinguish one flavor 
symmetry model from another through their different sensitivities to the radiative 
corrections. This observation makes sense in particular after the experimental errors 
associated with the neutrino mixing parameters are comparable with or smaller than 
the magnitudes of their respective RGE running effects. 

6. SEESAW-ENHANCED NEUTRINO DIPOLE MOMENTS 

The most popular mechanism of generating fi nite but tin y neutrino masses beyond 

190191 loo 9^ 94. 

the SM is the canonical seesaw mechanism 1 u | 11 J |1 , where the small neutrino 
masses are attributed to the existence of heavy degrees of freedom such as the right- 
handed Majorana neutrinos. In this elegant picture the 3x3 MNSP matrix U has 
a striking difference from the 3x3 CKM matrix V in the SM: it is not exactly 
unitary due to small mixing between light and heavy neutrinos as a result of the 
Yukawa interactions. The heavy neutrinos can be searched for at the LHC if the 
seesaw mechanism works at the TeV scale HS 7 | 138 | 139|14CH an( j |- ne un itarity of U 
can be tested in the future neutrino oscillation experiments. Another possibility 
is to look at the seesaw-induced non-unitary effects on the e lectromagnetic dipole 
moments (EMDMs) and radiative decays v i — » z/- + 7 11411. If the active neutri- 
nos acquire their respective masses in the seesaw mechanism, they should have the 
EMDMs through quantum loops. The fact that the Majorana neutrinos are their 
own antiparticles implies that they can only have the transition EMDMs between 
two different neutrino mass eigenstates in an electric or magnetic field. The rele- 
vant radiative decays of the heavier active n eutrinos wh ich may contribute to the 
cosmic infrared background in the Universe | 142 | 143 | 144 |^ are Q f p ar ticular interest 
in cosmology. 

6.1. Analytical Discussions 

We focus on the seesaw-induced non-unitary effects on the EMDMs and radiative 
decays of active Majorana neutrinos. The canonical seesaw mechanism is based on a 
simple extension of the SM in which three heavy right-handed neutrinos are added 
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and the lepton number is violated by their Majorana mass term 



20 21 22 23 24. 



-Cu = i l Y y HN K + \N-M R N R + h.c. , (52) 

where H = ia 2 H* with H being the SM Higgs doublet, £ L denotes the left- 
handed lepton doublet, N R stands for the column vector of three right-handed 
neutrinos, and M R is a symmetric Majorana mass matrix. After spontaneous 
SU(2) L ® U(1) Y — > U(l) em gauge symmetry breaking, H achieves its vacuum ex- 
pectation value (H) — v/y/2 with v ~ 246 GeV. Then the Yukawa-interaction term 
in C v yields the Dirac mass matrix M D = Y v v/\/2, but the Majorana mass term in 
C v keeps unchanged since right-handed neutrinos are the SU(2) L singlet and thus 
they are not subject to the electroweak symmetry breaking. The overall neutrino 
mass matrix turns out to be a symmetric 6x6 matrix and can be diagonalized 
through 

\K Mj \ M N ) 

where we have defined M v = Diag{m 1 , m 2 , m 3 } and M N = D'mg{M 1: M 2 , M 3 } with 
m i and M i being the physical masses of three light neutrinos v i and three heavy 
neutrinos N t (for i = 1,2, 3). The 6x6 unitary matrix U is decomposed as^J 

u=(ll)( A «)( x °). («) 



ZJ \S BJ \0 1 

where 1 denotes the 3x3 identity matrix, X and Z are the 3x3 unitary matrices, 
and A, B, R and S are the 3x3 matrices which characterize the correlation between 
the active or light neutrino sector (X) and the sterile or heavy neutrino sector (Z). A 
full parametrization of U in terms of 15 mixing angles and 15 CP-violating phases 
has been given before One may express the flavor eigenstates of three active 
neutrinos in terms of the mass eigenstates v i and N t . In the mass basis of charged 
leptons and neutrinos, the leptonic weak charged-current (cc) and neutral-current 
(nc) interactions read 

- C cc = JL jj— ^ {U mV ^ + R a . Nih ) W- + h.c] , 

+ [^(f/tii)..Ar. L + h.c] } Z^ , (55) 

where a runs over e, /i or t, U = AX is responsible for the flavor mixing of active 
neutrinos v i: and R measure s the strength of charged-current interactions of heavy 
neutrinos N i (for i = 1, 2, 3) 0321. A small deviation of U from X is actually char- 
acterized by nonvanishing R, as UU' = AA^ = 1 RR^ holds. The exact seesaw 
relation between the masses of light and heavy neutrinos is UM V U T + RM N R T = 0, 
which signifies the correlation between neutrino masses and flavor mixing parame- 
ters. 
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Fig. 4. The one-loop Feynman diagrams (and their charge-conjugate counterparts) contributing 
to the EMDMs of the Majorana neutrinos, where a = e, /i,t and i,j = 1, 2, 3. 




Fig. 5. The one-loop ■y-Z self-energy diagrams (and their charge-conjugate counterparts) associ- 
ated with the EMDMs of massive Majorana neutrinos, where / denotes all the SM fermions and 
i,j = 1,2,3. 

Let us consider the radiative v i — > + 7 transition, whose electromagnetic 
vertex can be written as 

rf.(0) = /iy (i (T^ft,) + e« K^^s) (56) 

for a real photon satisfying the on-shell conditions q 2 = and <z M £ p = 0. In Eq. (56) 
and /iy are the electric and magnetic transition dipole moments of Majorana 
neutrinos, and their sizes can be calculated via the proper vertex diagrams in Fig. 4 
(weak cc interactions). The j-Z self-energy diagrams in Fig. 5 (weak nc interactions) 
do not have any net contribution to and /Uy , but we find that they play a very 
crucial role in eliminating the infinities because the divergent terms originating from 
Fig. 4 are unable to automatically cancel out in the presence of the seesaw-induced 
non-unitary effects (i.e., R ^ and U 7^ X) unless those divergent terms originating 
from Fig. 5 are also taken into account. This observation is new. It implies that the 
non-unitary case under discussion is somewhat different from the unitary case (i.e., 
R = and WU = X^X = 1) discussed before in the literature E6 | 147 | 148 | 149 ] 
where the Feynman diagrams in Fig. 5 are forbidden and the divergent terms arising 
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from Fig. 4 can automatically cancel out. 
After a careful calculation, we arrive at 



ieG 



(m, i - m^) ^2 F « Re (^ai^a 



where 



2£ Q 



+ 



(57) 



(58) 



with £ Q = m^/M^ (for a — e, ^,r) denotes the one-loop fu nction Although this 
result is formally the same as that obtained in the literature | 146 | 147 | 148|149| they 
are intrinsically different as the seesaw-induced non-unitary effects on [i^ and 
were not considered in the previous works. To see how important such non-unitary 
effects may be, let us make two analytical approximations. First, F a ~ 3 (2 — £ a ) /4 
holds to a good degree of accuracy for £ Q <C 1. Second, U = AX ~ X — TX is also 
a good approximation for small non- unitary corrections to X, where ^B 



X 



c 12 c 13 

Coo Ci nS 



, 12 L 13 



'13 



'12 L 23 



12^13^23 °12°23 



'12 13°23 L 13' , 23 



s 12 s 23 C 12 S 13 C 23 



-12 s 23 s 12 s 13 c 23 c 13 c 23/ 



(l 



T = 



£*?* 



2 

fc=4 

6 

£ *lfc*2fe 

k=4 
6 



1 



\" „2 



2 ^ 

6 



"■2A- 



. a lfc s 3fe / , a 2fc a 3fc 
\fc=4 fe=4 



1 6 

2 ^ 

fc=4 



(59) 



with c,- 



and s 



(here ■ and are the mixing angles 



and CP-violating phases) . Note that the light-he avy n eutrino mixing angles 9 ik (for 
i = 1, 2, 3 and fc = 4, 5, 6) are at most of O(0. 1) USB sucn that the deviation of U 
from X is at the percent level or much smaller. Then we obtain 

: p a (u ai u* aj ) * -~ £ [(*)«, ( TX t 3 + ( Tx u w; 



£ 



|E[Wai( TX )a 
a 



(60) 



The first and second terms on the right-hand side of this equation correspond to the 
non-unitary and unitary contributions, respectively. While the former is suppressed 
by s 2 fc < 0(1O -2 ) (for i = 1,2,3 and k = 4, 5, 6) hidden in T, the l atter is suppressed 



by £ Q < 4.9 x 10" 4 (for a = e,/i, r) due to the GIM mechanism USB. We therefore 
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Fig. 6. Illustration of the seesaw-induced non-unitary effects on /i eff for three active neutrinos. 
The standard (unitary) results correspond to e uv = 0, and their uncertainties come from the three 
unknown CP- violating phases of X. 



draw a generic conclusion that the seesaw-induced non-unitary effects on e^- and a i 
can be comparable with or even larger than the standard (unitary) contributions. 
In this case the rates of radiative v i — > j/ ■ + 7 decays are given by 




with /i off = yj \fi i3 | 2 + \t i} | 2 for v i — » Uj + 7 being the effective EMDMs and 
/i B = e/(2m e ) being the Bohr magneton. The size of r„ + can be experi- 
mentally constrained by observing no emission of the photons from solar v e and 
reactor V e fluxes. More stringent constraints on u cS come from the Supernova 
1987A limit on neutrino decays and from the cosmological limit on distortions 
of the cosmic microwave background radiation (in particular, its infrared part): 
u eS < a few x 1CU 11 /i B | 152 | 153 | n ow more and more interest is being paid 
to the cosmic infra red background relevant to the radiative decays of massive neu- 
trinos HHH13II111 ; ^ i s desirable to evaluate /i cff and T v ^ u+ on a well-defined 
theoretical ground, such as the canonical seesaw mechanism under discussion. 



6.2. Numerical Illustration 

Figs. 6 and 7 illustrat e the numerical results of the non-unitary effects on fi eS and 
T u , respectively H3B Note that in this numerical analysis those small active- 
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Normal Hierarchy 
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Fig. 7. Illustration of the seesaw-induced non-unitary effects on r , for three active neu- 

trinos. The standard (unitary) results correspond to e uv = 0, and their uncertainties come from 
the three unknown CP- violating phases of X. 

steri le ne utrino mixing angles in Eq. (59) are constrained by present experimental 
dataQSQlas follows: 

T n < 5.5 x 1(T 3 , |T 21 | < 7.0 x 1(T 5 , 
T 22 < 5.0 x 10" 3 , |T 31 | < 1.6 x 10~ 2 , 
T 33 < 5.0 x 10~ 3 , |T 32 | < 1.0 x 10~ 2 . (62) 

All s ik in T (for i — 1,2,3 and k = 4,5,6) are positive or vanishing. The CP 
phases S ik are all allowed to vary from zero to 2tt, but they must satisfy the above 
constraints together with the relations UW+RK* = 1 and UM V U T +RM N R T = 0. 
To assure that radiative corrections to the masses of three light neutrinos (via the 
one-loop self-energy diagrams involving the heavy neutrinos) are sufficiently small 
(e.g., smaller than 0.5 meV) and stabl e, we sim ply assume that the masses of three 
heavy neutrinos are nearly degenerate | 154 | l5 5] an j no ^ morc than 0(1) TeV. This 
assumption implies that the results shown here are for a limited and safe parameter 
space of the TeV seesaw mechanism, but it is instructive enough to reveal the salient 
features of the non- unitary effects on the effective EMDMs MeffC^j — ^ v j + 7) ano - 
the radiative decay rates T u +7 . 

To present our numerical results in a convenient way, let us define 

~\ 1/2 

(63) 



6 

J2 (4 k + 4k + 4k) 

,fc=4 



which measures the overall strength of the unitarity violation of U, and e uv € 
[0,0.15) is reasonably taken in our calculations. Namely, we allow each s ik (for 
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i = 1,2,3 and k — 4,5,6) to vary in the range < s ik < 0.15. The numerical 
dependence of /^(z^ — > Vj + 7) and r„ +7 on e uv is shown in Figs. 6 and 7, 
respectively. Some discussions are in order. 

(1) Switching off the non-unitary effects (i.e., e uv = 0), we obtain the effective 
electromagnetic dipole moments 

r (0.8 - 3.0) x 10~ 25 fi B (i> 2 v 1 + 7) , 
Mcff - { (0-8 ~ 1.5) x 10~ 24 Mb (1/3 -> ^ + 7) , (64) 
[(1.1^2.1)x10- 24 Mb K^^2+7), 

for the normal mass hierarchy with m 1 ~ 5 meV; and 

r (0.01 ~ 2.0) x 10~ 24 Mb (f a ->• fi + 7) , 
Mcff ~ ^ (0.8 ~ 1.5) x 10- 24 Mb (1/3 -> v x + 7) , (65) 
[ (1.3 ~ 2.0) X 10" 24 /z B 

for the inverted mass hierarchy with to 3 ~ 5 meV, where the uncertainties mainly 
come from the unknown CP phases <5 12 , <5 13 and <5 23 . Such standard (unita ry) result s 
are far below the observational upper bound on fi cS (< a few xl0~ n fi B I I I), 
but they serve as a good reference to the non-unitary effects on /i cff being explored. 

(2) Figs. 6 and 7 clearly show that /j, eS and Y v ^ v + ^ can be maximally enhanced 
by a factor of O(10 2 ) and a factor of O(10 4 ), respectively, in particular when e uv 
approaches its upper limit as set by current experimental data. The magnitude 
of At ef j(^2 ~ ^ ^1+7) m ay be strongly suppressed in the inverted neutrino mass 
hierarchy. The reason is rather simple: on the one hand, m 1 ~ m 2 holds in this 
case, and thus e 12 oc (m 2 — must be very small; on the other hand, fi 12 depends 
on lm(U al U* 2 ), so it can also be very small when the CP-violating phases are 
around zero or tt. This two-fold suppression becomes severer for the decay rate 
r„ _y v j^, because it is proportional to (m 2 — Wi) 3 /^ (^2 — > V\ + 7). 

(3) The results of /x eff and r„ +7 are sensitive to the absolute neutrino mass 
scale for both normal and inverted mass hierarchies. For instance, H c g(v 2 ^1+7) 
and /i eff (f 3 — > 17 + 7) get enhanced when m 1 changes from zero to 5 meV in the 
normal mass hierarchy; while fj, cS (v 1 — > v 3 + 7) and /x eff (i/ 2 — > ^ 3 +7) are enhanced 
when m 3 changes from zero to 5 meV in the inverted mass hierarchy. This kind of 
sensitivity is not so obvious if one only takes a look at the expressions of ■ and 

in Eq. (57). The main reason is that a change of m 1 or m 3 requires some fine-tuning 
of the active-sterile neutrino mixing angles and CP-violating phases as dictated by 
the exact seesaw relation UM U U T + RM N R T — 0, leading to a possibly significant 
change of /i eff . The dependence of T v ^ v +1 on the absolute neutrino mass scale is 
somewhat more complicated, as one can see from Eq. (61). 

(4) The CP-violating phases play a very important role in fitting both the exact 
seesaw relation and Eq. (62). If the heavy neutrino masses M i are not suppressed, 
then an appreciable value of e uv requires some fine cancellations in the matrix 
product RM N R T such that sufficiently small m i can be obtained from UM^U 7 = 
—RM N R T . On the other hand, we remark that it is actually unnecessary to require 
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M i to be around or above the electroweak scale. The seesaw-induced non-unitary 
effects on fj, eS and r„ __ j . v + can be significant even if one allows one, two or three 
heavy neutrinos to be relatively light (e.g., at the keV mass scale). Such sterile 
neutrinos are interesting in particle physics and cosmology. Note that it is easier to 
satisfy the exact seesaw relation with an appreciable value of e uv by arranging M i to 
lie in the keV, MeV or GeV range. This kind of low-scale seesaw scenarios 1^ might 
be technically natural, but they have more or less lost the seesaw spirit. Of course, 
sufficiently large M i and sufficiently small 9 ik can always coexist to make the seesaw 
mechanism work in a natural way, but in this traditional case the non-unitary effects 
are too small to have any measurable consequences at low energies. 

It is also worth pointing out that the seesaw-induced non-unitary effects on fj,^ 
and are rather different from the case of making a naive assumption of the flavor 
mixing between three active neutrinos and a few light sterile neutrinos ^Sl. The 
latter can directl y bre ak the unitarity of the 3x3 MNSP matrix U and then lift the 
GIM suppression USD associated with fj,^ and . This kind of non- unitary effects 
are not constrained by the seesaw relation, and thus they are more arbitrary and 
less motivated from the point of view of model building. 

The effective electromagnetic dipole moments of three neutrinos and the rates of 
their radiative decays can be maximally enhanced by a factor of O(10 2 ) and a factor 
of O(10 4 ), respectively, no matter whether the seesaw scale is around or below the 
TeV energy scale. This observation is new and nontrivial, and it reveals an intrinsic 
and presumably important correlation between the electromagnetic properties of 
neutrinos and the origin of their masses. Such a correlation may even serve as a 
sensitive touch-stone for the highly-regarded seesaw mechanism. 

7. SUMMARY 

After the Daya Bay measurement of the smallest neutrino mixing angle 8 13 , it is 
natural to ask where we are standing and where we are expecting to go in neutrino 
physics. We have tried to answer these two questions from a phenomenological point 
of view in this review paper, although our answers are incomplete and full of con- 
jectures. To be specific, we have given a fast overview of some fundamental neutrino 
properties and paid particular interest to the flavor issues of charged leptons and 
neutrinos, including the mass spectrum, flavor mixing pattern and CP violation. 
We have gone into details of possible lepton flavor structures by describing two 
useful phenomenological strategies and giving a number of typical examples. The 
impact of large 9 13 on the running behaviors of other flavor mixing parameters has 
been discussed in the framework of the MSSM. We have also illustrated the seesaw- 
enhanced electromagnetic dipole moments of three Majorana neutrinos based on a 
viable TeV seesaw scenario. 

If only the SM particles are taken into account and the massive neutrinos are 
assumed to be the Majorana particles, we are then left with 29 fundamental pa- 
rameters in Nature, as described by the so-called Fritzsch-Xing plot in Fig. 8. The 
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Fig. 8. The Fritzsch-Xing plot for 29 fundamental parameters in Nature, describing four kinds 
of interactions (G N , a, a w and ce s ), six quark masses (m u , m c , m t and m d , m s , m b ), six lepton 
masses (m e , m„, m T and ra 1 , m 2 , ni 3 ), three quark flavor mixing angles (l? 12 , #13 and $23)' three 
lepton flavor mixing angles (0 12 , ^3 and $23)' two CP- violating phases in the quark sector (S q 
and 8), three CP- violating phases in the lepton sector (<5 ; , p and <r), the Higgs mass M H and the 
H/-boson mass M w . 



determination of 6* 13 and M H in 2012 is a milestone in particle physics. But the 
effective strong CP-violating phase 5 in the quark sector and the three weak CP- 
violating phases (i.e., <5 ; , p and a) in the lepton sector remain unknown. Moreover, 
the absolute mass scale of three neutrinos and their mass ordering have not been 
determined. We hope that various precision neutrino experiments can help pin down 
the relevant parameters in the lepton sector and then shed light on the flavor dy- 
namics in the foreseeable future. 

Of course, the picture in Fig. 8 may be too simple and too naive because we are 
not sure whether some of those "fundamental" parameters are really fundamental 
or not. New degrees of freedom, such as the sterile neutrinos or the supersymmetric 
particles, might be discovered and make the flavor sector much messier. If the under- 
lying flavor theory is regarded as an animal, one has no idea whether it is a donkey 
or an elephant or something else. In this sense we are blind today and have to make 
a lot of experimental and theoretical efforts to identify its nose, eyes, ears, legs and 
so on in order to make sure what animal it is. The road behind has repeatedly told 
us that the road ahead is always challenging, but it is always exciting. 



This review paper is essentially based on the plenary talk given by one of us 



November 28, 2012 3:16 WSPC/INSTRUCTION FILE SUSY 



Theoretical Overview on the Flavor Issues of Massive Neutrinos 33 

(Z.Z.X.) at the SUSY 2012 conference. The work of S.L. is supported in part by 
the National Basic Research Program (973 Program) of China under Grant No. 
2009CB824800, the National Natural Science Foundation of China under Grant 
No. 11105113 and the Fujian Provincial Natural Science Foundation under Grant 
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